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Abstract 

We consider the Farey fraction spin chain, a one-dimensional model defined 
on (the matrices generating) the Farey fractions. We extend previous work on 
the thermodynamics of this model by introducing an external field h. From 
rigorous and more heuristic arguments, we determine the phase diagram and 
phase transition behavior of the extended model. Our results are fully consistent 
with scaling theory (for the case when a "marginal" field is present) despite the 
unusual nature of the transition for h — 0. 

Keywords: phase transition, Farey fractions, spin chain, 



'Electronic address: jan.fiala(at)umit. maine.edu, kleban(at)maine. edu; fax: (207)-581-2255; phone: (207)-581-2258 



1 



I. INTRODUCTION 



Phase transitions in one-dimensional systems are unusual, essentially be- 
cause, as long as the interactions are of finite range and strength, any pu- 
tative ordered state at finite temperature will be disrupted by thermally 
induced defects, and a defect in one dimension is very effective at destroying 
order. Despite this, there are many examples of one-dimensional systems 
that do exhibit a phase transition. The Farey Fraction Spin Chain (FFSC) 
is one such casejwhich has attracted interest from both physicists and 
mathematicians 

@,H0. ( Since this work 

uses some methods that may be 
unfamiliar to the latter, we include a paragraph at the end of this section 
outlining our results from a mathematical viewpoint.) 

One can define the FFSC as a periodic chain of sites with two possible 
spin states (A or B) at each site. This model is rigorously known to exhibit 
a single phase transition at temperature (3 C = 2 |]|. The phase transition 
itself is most unusual. The low temperature state is completely ordered 
. In the limit of a long chain, for j3 > /3 C , the system is either all A 
or all B. Therefore the free energy / is constant and the magnetization m 
(defined via the difference in the number of spins in state A vs. those in 
state B) is completely saturated over this entire temperature range. Thus, 
even though the system has a phase transition at finite temperature, there 
are no thermal effects at all in the ordered state. The same thermodynamics 
occurs in the Knauf spin chain (KSC) 00, ED, EH, to which the FFSC is 
closely related. 

At temperatures above the phase transition (for f3 < /3 C ), fluctuations 
occur, and / decreases with (3. Here the system is paramagnetic, since 
(when the external field vanishes, see below) there is no symmetry-breaking 
field. Thus as the temperature increases m jumps from its saturated value 
in the ordered phase to zero in the high-temperature phase 0, Q (see 
Fig. GJ. (The KSC behaves similarly.) 

One-dimensional models with long-range ferromagnetic interactions 
are known to exhibit a discontinuity in m at /3 C , but in these cases the jump 
in m is less than the saturation value. 

The discontinuity in m might suggest a first-order phase transition, but 
in our model the behavior with temperature is different. In previous work, 
we proved that as a function of temperature, / exhibits a second-order 
transition, and the same transition occurs in the KSC and the "Farey tree" 
multifractal model 0). 

In beginning the research reported here, our motivation was to see 
whether the phase transition in the FFSC, which seems to mix first- and 
second-order behavior, is consistent with scaling theory Indeed, as will 
be made clear, it is, in the "borderline" case when a marginal variable is 
present. In order to see this, we extend the definition of the FFSC to include 
a finite external field h. We then determine the phase diagram and free 
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FIG. 1: Free energy and magnetization vs reduced temperature t = — 1 



energy as a function of (3 and h, using both rigorous and renormalization 
group (RG) analysis. 

In the following, section [Q] defines the model. Then, in section II I II we 
prove the existence of the free energy / with an external field, and evaluate 
/ for temperatures below the phase transition. In section IIVI we employ 
renormalization group arguments to find the free energy and phase dia- 
gram for temperatures above the phase transition. Section IVI considers a 
simple model that has very similar thermodynamics but is completely solv- 
able. Section IVll summarizes our results. In the Appendix we present some 
arguments needed to prove the existence of /(/?, h) in section ITTT1 

Since our results may be of interest to mathematicians who are unfa- 
miliar with some of the physics employed herein, we pause to include a 
description of them from a more mathematical point of view. Section [H] 
defines the model and the quantities of interest. More specifically, the par- 
tition function Zn is a two-parameter weighted sum over the (matrices 
defining the) Farey fractions, and the free energy / then follows from the 
limiting procedure defined in Q . The main goal of our work is to find the 
analytic behavior of / as a function of the real parameters (3, the inverse 
temperature (so f3 > is implicit), and h, the external field. Regions of pa- 
rameter space for which / is analytic are (thermodynamic) phases, and the 
lines of singularities that separate them are phase boundaries. In section 
lllll we prove that h) exists, and compute it exactly at low temperature 
(for (3 > f3 c ), which constitutes part of the ordered phase. Section llVl uses 
renormalization group methods to determine / at high temperatures (for (3 
near (3 C and (3 < (3 C ). Since this method is not rigorous, from a mathemat- 
ical point of view the results should be regarded as conjectures. The main 
conclusions are the form of the free energy in the high-temperature phase 
1)291 I30fl , the equation for the phase boundary ()31l 132)1 and the change in 
magnetization m — —df/dh 1)340 and entropy s = (3 2 df/d/3 across the 
phase boundary. We also find that the ordered phase, with / = ex- 
tends to (3 < (3 C when h is sufficiently large (see Fig. (3J . Section IIV CI 
gives predictions for the behavior of Zn as N — > oo near the second-order 
point {(3 = (3 C and h — 0). This is related to some work in number theory, 
but unfortunately not yet directly. Section examines an exactly solvable 
model with certain similarities to the FFSC. 
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II. DEFINITION OF THE MODEL 



The FFSC consists of a periodic chain of N sites with two possible spin 
states (A or B) at each site. The interactions are long-range, which allows 
a phase transition to exist in this one-dimensional system. Let the matrices 



M N := JjA 1 -"^', 0*6 {0,1}, (1) 

where A := and B :— ( 01 ) and the dependence of Mn on {cr;} has 

been suppressed. The energy of a particular configuration with N spins in 
an external field h is given as 

E N := ln(Tjv) + h ( 2^~Vi - N] with T N := Tr(Mjv). (2) 



Thus our partition function is 

Z N (J3,h) = Tr(M N )-^e-^ h ( 2 ^=^'- N ). (3) 
{*<} 

This definition extends the Farey fraction spin chain model to non- 
vanishing external field h. Given the nature of the low-temperature h = Q 
system, it is natural to introduce h in this way. 
The free energy is defined as 

/(A h) := ^ lim ^E^R. (4 ) 

The existence of the free energy f{(3, h) follows from simple bounds using 
/(/3,0) (see section Hill below) . 

The definition of the FFSC is somewhat unusual. The partition function 
is given in terms of the energy of each possible configuration, rather than 
via a Hamiltonian. In fact, there is no known way to express the energy 
exactly in terms of the spin variables a . Further, numerical results indicate 
that when one does, the Hamiltonian has all possible even interactions (and 
they are all ferromagnetic), so an explicit Hamiltonian representation, even 
if one could find it, would be exceedingly complicated. 

Note that for h = there are two ground states with energy E = In 2. 
The other 2^ — 2 states have energy In N < E < Nc, where c is a constant. 
Therefore the difference between the lowest excited state energy and the 
ground state energy diverges as N — > oo. 

The phase transition in this system Q occurs in the following way. Divide 
the partition function into two terms, one due to the two ground states, 
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and the other (call it Z'), due to the remaining 2 N — 2 states. The system 
remains in the ground states, and Z' — > as N — > oo, until the temperature 
is high enough that Z' diverges with N. In section^we examine a simple 
model that also exhibits this feature, but is completely solvable. 

Our results also apply to the KSC. which has the same thermodynamics 
as the FFSC model at h = (see 2J). An external field may be included 
in the KSC in exactly the same way as described above for the FFSC. The 
"Farey tree" model of Feigenbaum et. al. ^3 a l so nas the same free energy, 
but it is not clear how to incorporate a field h. Our finite-size results (see 
section flV CI) do apply when h — 0, however. 



III. FREE ENERGY WITH AN EXTERNAL FIELD 

In this section we show rigorously that f(/3, h) exists and that 

f(J3,h) = -\h\, (5) 

for ft > ft c . 

For h > it is easy to see (from Q) that 

Using the definition of the free energy then gives 

-h>f(ft,h)>f(ft,0)-h, (7) 

where f(ft, h) is understood to be defined via ((2J|. Now f(ft, 0) is rigorousl 
known to exist In addition, we know that f(/3, 0) = for (5 > {3, 
which implies for h > (h < follows similarly). 

To see that /(/?, h) exists for the range < (5 < f3 c we proceed as 
follows (actually, our argument applies for all (3 > 0). We first show that 

log Z N + 1 _ log Z A 
JV+1 N 

Now 



as N — > oo. The result then follows by use of JfJJl. 



N log Zjv + i — N log Zjy — log Zjm 



N(N + 1) 



< 



log Z N+ i/Z 



N 



N + l 



1 



logZjv 



N 



and we see by © and the existence of /(/3, 0) that the second term 



l 

iV+l 



log -Z^ 



N 



< for some finite constant K. In the appendix we show 

that 2~Pe~P <{h \ < g +1 < 2e^ h \ which completes our proof of the existence 
of the free energy for all ft > and h G M. 



We also know rigorously that f(t, 0) 



'int 



., where c > 0, t 



£§■ — 1, for t > (see Fig. [TJ. It follows that f(t,h) must have at least 
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one singularity between the regions with low and high temperatures, i.e. a 
phase transition from the ordered to the high-temperature phase. 

Since we can not calculate f(j3,h) exactly for j3 < f3 c (except for h = 
and [3 — > (3 C ), we use another method, in the next section, to examine the 
thermodynamics. 

IV. RENORMALIZATION GROUP ANALYSIS 
A. Mean field theory 

In mean field theory one assumes that there is an expansion of the free 
energy of the form 



where M is the magnetization and the "constants" a, b, u and g are weakly 
dependent on the reduced temperature t (defined at the end of section 
llllll and external field h. Note that u > is required for stability, and 
b > 0, g > in the high-temperature phase. (The possibility that g = is 
ruled out below.) 

Minimizing (JSJ with respect to M, one obtains the free energy and mag- 
netization in mean field approximation. Explicitly 

1. for t > and h ^ the magnetization 



(note the limiting cases Mo ~ for h = and Mo ~ 1/6((?/i/m) 1 / 3 
for t = 0) 

2. for t < and h ^ 0, but h sufficiently small, the magnetization 



(however when + 4 (||) > 0, M is given by the t > for- 

mula). We include this second case only for completeness. Since our 
system is completely saturated at low temperatures this result is not 
employed in our analysis. 

In the following we use the first result in an RG analysis. 



f MF = a + MM 2 + uM 4 - ghM + ... 



(8) 
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B. Renormalization group analysis 



We assume two relevant fields (t and h) and one marginal field (u). 
These assumptions are reasonable, since our model has an Ising-like ordered 
state, the interactions are (apparently) all ferromagnetic, and there is a 
logarithmic term in the free energy. 

The infinitesimal renormalization group transformation for the singular 
part of the free energy is 



f s (t,h,u) = e' M f s (t(£),h(£),u(£)). 



(9) 



Because of the marginal field u, the analysis is somewhat more complicated 
than otherwise. We follow the treatment of Cardy (H3i see also Wegner 
14] ). The RG equations take the form 



du/d£ — -in 2 + . . . 
dt/dl — y t t — z t ut + . . 
dh/dl = yhh — Zhuh + 



(10) 

(11) 
(12) 



where we keep only the most important terms. The omitted terms are 
either higher order or go to zero more rapidly with I than those included. 
From (Unj we find (note t = t(0), h = h(0), u = u(0)) 



nil) = U{0) 

{> 1 + 314(0)*' 

Both t and h have the same functional form, namely 

Zt 



and 



ln(t(4,)/t(0)) = ytlo - - ln[l + xu(0)e o ] 



\n(h(e o )/h(0)) = y h t - ^ ln[l + xu(0)£ ] 

x 

where £q is such that t(£o) = 0(1) or h(£o) — 0(1). From l|14|l we can write 

to' 



(13) 



(14) 



(15) 



yt t xy t 



1 



— u In 

Vt t 



(16) 



where we assume to/t 1. This result together with @ gives us 



fa(t,h, u) 



t 


vt 







1 



X t 

— u In — 

yt t 



f s (t(to),h(£ ),u(£ )). (17) 



Since the free energy on the rhs is evaluated at ^o, which is far from the 
critical point, it can be calculated from mean field theory. Above the critical 
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temperature [t > 0) with small external field h (h(£o) <C t(£o)) we obtain 
for the free energy 



Hgh(e )) 2 
mt(£ ) 



(18) 



The relation between h(£o) and t(£o) follows from (|14|l and l|15|) . Eliminat- 
ing /i(^o) allows us to rewrite (fTH)l as 



/« ~ a - 



1 , a; t 
1 H u In — 

yt t 



2y h 



35 2 



16bt(4) 



Substituting the result into (0 gives two terms, 



t 


vt 


to 





1 , x 1 to 

1 H u In — 



and 



X t 

— tun — 

t 



Vt 



3.9 2 



166t(4) 



(19) 



(20) 



(21) 



The first term can be compared with the exact result at h = (see 
section [J. It follows that 



- 1 - £i 



(22) 



The second term gives us the dependence on external field. Eliminating 
t(£o) instead of h(£o) we obtain 



1 




d ,vt 
Vh Vh 


t 


h 





1 I x 1 to 

1 H u m — 

2/t * . 



3(gMl )) s 
166 



(23) 



Equating the two expressions l|21|) and (|23() for the same term in the free 
energy gives us the RG eigenvalues 



d _ d 

yt yh 



1. 



(24) 



where d is the dimensionality of the system. This is of course one for 
our model, but since none of our results require setting d = 1 we leave it 
unspecified. 

Finally we can write down the singular part of the free energy for the 
high-temperature phase 



Since / < for h = in this phase, (|25|l implies that a < 0. 

For the ordered phase we know rigorously that the free energy has no 
temperature dependence for h = 0. The spins are all up or all down. When 
we add an external field it will break the symmetry and all the spins will 
be oriented in the field direction. Thus the free energy at £o is 

f s (t(e ),h(e Q ),u(e )) = -\h\(e ). (26) 

Proceeding as in the derivation of (|17(l from and Ijl4(l we get 



/«(*, h,u) 



h 


Vh 


h 





x ho 

1 H urn — 

Vh h _ 



fs(t(e ),h(£ )Mto)), (27) 



using JIHJl, (J2EI and (jSJl then give 



/ s (t,M)~-N 



1 H u In — 

2//!. h 



(28) 



Because the magnetization in the ordered state is completely saturated the 
logarithmic correction must vanish. Therefore Zh = 0. 

Thus the asymptotic form for the free energy of the high-temperature 
state is 



fs(t,h, u) 



X t Q 

— urn — 

Vt t 



t 



X t 

— u In — 

Vt t 
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We can recast this result more suggestively as 

f s (t,h,u) ~f s (t,0,u)- h 2 X (t,0,u), 



(29) 



(30) 



where x — —d 2 f/dh 2 is the susceptibility. Note that \ ~ V/s which 
is consistent with scaling theory, since (using (J23J), / ~ t 2 ~ a — t d l yt = 
t while x ~ ^ 7 = t^y^/y* = t~ l . This relation holds regardless of 
whether we set the dimensionality d — 1 or not. In addition, the coefficient 
of j^t for the free energy at h = and t — > 0, i > is known exactly 
[TUnH' so that the combination of constants t ^j° M may be determined. 

The phase boundary is given by the continuity of the free energy. Now 
we expect the ordered phase to exist for (3 < (3 C if h is large enough (this is 
reflected in the assumption of two relevant fields-if another phase intervened 
there would be more). Thus one must equate the two expressions for /. One 
finds that the phase boundary between the ordered and high-temperature 
phase, close to the critical point, follows 



1*1 



t 



'hit /t '' 



(31) 
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I . Since / is quadratic in h in the high- 
temperature phase, there are in general two solutions with h > 0. However, 
the one at larger h is not physical since it gives rise to a magnetization 
m > 1 and violates the convexity of the free energy as well, so we employ 
the other. 



where k 



Zxug 2 



: »"/- 
ibt„ 



FIG. 2: Phase diagram 



In order to find the change in magnetization across the phase boundary 
we use (13 111 with constants included 



\h\ 



In* 



8by t 
3xug 2 



3ag 2 
4Mq 



(32) 



In arriving at l|32|) , we (as mentioned) chose the root that makes m < 1 in 

o 2 

the high-temperature phase. Note that in the limiting case that = — 1, 
?7i = 1 but the two roots coincide. 
Now from (jSHJ 



X t 

— u In — 

Vt t _ 



(33) 



Eliminating the external field using (|32|l . and since the magnetization in 
the ordered phase takes the values m ~ ±1, we find 



Am 




(34) 
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Note that to is a constant of order one and recall that a < 0, thus on the 
phase boundary the discontinuity in magnetization is constant (and less 
than one), at least close to the second-order point (we argue below that 
g = is not possible in this model). Now we can look at the change in 
entropy (per site) s = f3 2 df/d/3 across the phase boundary. We get 



x t 
— u In — 

Vt t 



a 4b 



1 - 4 1 



3ng 2 
Abto 



(35) 



These results show that the phase transition is first-order everywhere except 
at h = 0. 

3a 2 

In the limiting case when -gjrj— = — 1, already mentioned, one finds that 
both Am = and As = 0. However, it is easy to see that both the 
susceptibility \ an d the specific heat will have a discontinuity across the 
phase boundary. 

Note that the magnetization change given by l|34|l exhibits a kind of "dis- 
continuity of the discontinuity" , in that its limiting value as one approaches 
the second-order point is not the same as its value at that point. This is 
not the case for the entropy change, or for these quantities in the model 
examined in section Ivl 

Finally, we argue that g = is not possible in the high-temperature 
phase. Since the second derivative of / with respect to h at h — is pro- 
portional to both g and the susceptibility Xi h suffices to demonstrate that 
\ > 0. It is straightforward to show that x is proportional to S^ =1 (siSj) 
where the spin variables Sj := 2(7j — 1, Si G { — 1,1} (cf. (JTJ), and the 
angular brackets denote a thermal average. Now the j ' = 1 term in this 
sum is 1, and due to the ferromagnetic interactions in the spin chain, the 
remaining terms cannot be negative. Note that this argument is not com- 
pletely rigorous, since for the FFSC we only have numerical evidence that 
the interactions are all ferromagnetic. The KSC, on the other hand, is 
known to have all interactions ferromagnetic 0, so that (s\Sj) > follows 
from the GKS inequalities. 



C. Finite-size scaling 



We can use our results to make some predictions about finite-size (i.e. N 
large but N < oo) effects on the thermodynamics. We make the standard 
assumption that the size of our spin chain is a relevant field with eigenvalue 
1. Of course, since our system has long-range interactions the validity of 
finite-size scaling may be questioned [13|. but it is still interesting to see 
the results. The treatment is the same as in the case of the relevant fields t 
and h. The renormalization equation for the inverse size I := N^ 1 is then 

dl/dl = l-ziul + ... . (36) 
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Thus we get 



No 


d r 


N 





1 + x u In ■ 



§-] * A(t(4),M4)X4),A^Vo)). 

(37) 

Note that we do not know the ratio Zi/x, however (|37() gives the form we 
should observe. More succinctly, for large N, this result predicts that for 
small t and h 



(38) 



There is related work in number theory by Kanemitsu 01 ( c f- a l so 
[Tsj)- This paper studies moments of neighboring Farey fraction differ- 
ences, which are similar to the "Farey tree" partition function At 
h = 0, the latter has the same thermodynamics as the FFSC 0. However, 
[l7j uses a definition of the Farey fractions that, at each level, gives a subset 
of the Farey fractions employed here, and none of the moments considered 
corresponds to f3 = 2 (the point of phase transition). It is interesting that, 
despite these differences, terms logarithmic in N appear. More specifically, 
the sum of mth (integral) moments of the differences goes as 



(In AO 



AT" 



O 



(In N) 



h(rri) 



(39) 



for m > 2, with g(2) — 1, <j>(3) = 2 and g(m) = 3 for m > 4, and h{m) = 1 
for 2 < m < 4, /i(m) = for m > 5. Now if all the Farey fractions 
were included (1391) would apply to the Farey tree partition function with 
/ 3 = 2m(cf. UllJ) so that m > 2 would correspond to (3 > 4. It would be 
interesting to determine whether l|39|l applies to the Farey tree partition 
function despite this difference, or to extend l|39() to m = 1 to see if it is 
consistent with l|38|) . 



V. 1-D KDP MODEL WITH NONZERO EXTERNAL FIELD 

In this section we consider the one-dimensional KDP (Potassium dihy- 
drogen phosphate) model introduced by Nagle This model's thermo- 
dynamics and energy level structure are similar to the Farey fraction spin 
chain, but it is easily solvable. Comparison of the two models thus sheds 
some light on the FFSC. 

The KDP model exhibits first-order phase transitions only. The origin 
of the phase transition is infinite rather than long-range interactions. 

The one-dimensional geometry of the model is illustrated in Fig. |3J It 
consists of N cells, and each cell contains two dots. Each dot represents a 
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A 



B 



C 



D 



FIG. 3: KDP 



proton in a hydrogen bond in the KDP molecule. Dots can be on the left or 
the right side of a cell. The energy of a neighboring pair of cells depends on 
the arrangement of dots at their common boundary. Only configurations 
with exactly two dots at each boundary (e.g. A, B and D in Fig. I3J) are 
allowed, any other configuration (e.g. C in Fig. |3J) has (positively) infinite 
energy and is therefore omitted. Of the allowed configurations, only two 
energies occur, (when there are two dots on the same side of a boundary, 
as in Fig.[3]D) or e (when the dots are on opposite sides, as in Fig. [31 A or 



Let there be N cells in a chain with periodic boundary conditions. Then 
there are two kinds of configurations with finite energy. In the first type 
of configuration, each cell has two dots on the same side. There are two 
such configurations and the total energy of each is 0. In the second type of 
configuration, each cell has one dot on the left and one on the right. There 
are 2 N such configurations and the total energy of each is Ne. Thus, the 
partition function is simply 



It follows immediately that / = for fle > In 2 and / = e- ^ for (3e < In 2. 
Thus the temperature of the (first-order) phase transition is T c = e/(ln2) 
and there is a latent heat with entropy change As = In 2. Clearly, the 
phase transition mechanism is a simple entropy-energy balance. At low 
temperatures, the ground state energy gives the minimal free energy, while 
in the high-temperature phase the extra entropy of the additional states 
gives a lower free energy. 

Next, define the magnetization m as the number of sides of cells with 
both dots on one side divided by the number of cells N. Then m = 1 for 
(3 > (3 C and m = for [3 < (3 C (so that Am = 1 at the phase transition), 
just as in the FFSC model. 

Following the above definition of the magnetization, we introduce an 
external field h by adding an energy ±h/2 to each dot, according to whether 
it is on the right or left side of the cell. This gives the extra energy of an 
external field acting along the chain. Then the new partition function has 



B). 



2 + 2 N exp(-/3Ne). 



(40) 
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the form 



Z N (f3, ft) = 2 cosh(/?iV/i) + 2 N exp(-/3iVe). (41) 

In the ordered phase Z — > exp(±/3iVft). Thus, the free energy / = qp/i, 
where the plus sign is for ft, > and the minus sign for ft < 0, exactly as in 

and we 



the FFSC. For the high-temperature phase Z — » exp 7V(ln2 — /3e) 

lnj 

/3 



get the same free energy as when ft = 0, / = e — The phase boundary 



is given by ft. = ±et (see Fig. 0J, where i = % — 1 as before. Note the 
resemblance to the FFSC phase diagram (Fig. [21). Here, as (3e — > In 2, ft — > 
as it should, while for /3 — > the field ft. — ► ln2//3. The entropy per site 
vanishes everywhere in the ordered phase, while for the high-temperature 
phase s = In 2. Thus, this model has a non-zero latent heat and the phase 
transition is first-order everywhere. Note that the change in magnetization 
is Am = 1 everywhere along the phase boundary between the ordered state 
and the high-temperature state. 



FIG. 4: Phase diagram 

Now for ft = 0, the FFSC has two ground states with all spins up or all 
spins down and energy independent of length N, just as in the KDP model. 
Then, in addition, the FFSC has 2 N — 2 states with energies between IniV 
and Nc, for some constant c. On the other hand, the KDP model has just 
one energy (Ne) for the 2 N states corresponding to the 2^ — 2 states of the 
Farey model. This might suggest that the states with energies close to In N 
are responsible for the logarithmic factor in the Farey free energy, and thus 
shift the phase transition from first to second-order (for ft = 0). For ft ^ 
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the energy of the In N states is shifted by the field h to order N, and the 
phase transition becomes first-order. However the mechanism of the FFSC 
phase transition may be more subtle. The "density of states" (number of 
configurations with a given energy) for the FFSC not well-behaved. In fact 
it is known rigorously that this quantity, summed over all chain lengths, 
has a limit distribution [3j. 

Note that the free energy just derived is independent of h in the high- 
temperature phase. Since this is not what we found for the FFSC, we 
consider another way to introduce an external field h into the KDP model. 



field 



1 

Si : + 1 
ft: 



2 
-1 





3 


+ 1 



4 

-1 



FIG. 5: notation 



As before we have four different states for each cell. We index them with 
spin-one variables U and Si (sj, U £ {0, +1, —1} ) in each cell as in Fig. El 
Then the energy (for h = 0) can be written 



JV-l 



H - e ^2 



(42) 



»=i 



(assuming, in the sum, that the infinite energy contributions are omitted). 
The conditions s, + U = ±1 and Sjij = define the allowed states. We 
define the magnetization per site as 



1 N 



(43) 



Note that this definition gives a positive (negative) contribution if the upper 
dot in a given cell is on the right (left). (Note also that m 2 



N 2 



Ei=iOi+ 

ti) 2 + -faEi*M+ti)(«j+tj) = &Ejl 2 («i+*i)( fl j + *j) + VW-) Hence 
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we can include an external field as follows 



H = H -h^2(s t + U) = H - hNm. (44) 

i=l 

Thus 

Z{(3, h) = e 0Nh + e-? Nh + e~^ N [2 cosh(/3/ l )] Ar , (45) 
and the free energy in high-temperature phase becomes 

Mh)=e- H2CO f m) (46) 

or for small h 

f~-te- n ! U2 , h\ (47) 

J 2e(i + l) v ' 

with t = — 1 as above. The phase boundary is given by 



ph = In (2 cosh(/3/i)) - (3e. (48) 



For |3/i<Cl and /i > 0, using /3 C = — , this gives 

e In 2 t „ . o , , 
ft = et+ ~^-t 2 + 0{t 3 ), (49) 

The phase diagram near the critical point is very close to the previous one 
(see Fig.0J. The magnetization in the ordered phase is again independent 
of temperature, i.e. m = ±1. In the high-temperature phase we have 
m = tanh(/3/i). Thus the magnetization change across the phase boundary 
close to the critical point is Am = 1 — <ln2. The transition is again first- 
order, with the entropy change As = In 2(1 — ^i 2 ). Results for h < 
follow immediately by symmetry. 



VI. SUMMARY AND COMMENTS 



In this paper, we have extended the definition of the Farey fraction spin 
chain to include an external field h. From rigorous and more heuristic 
arguments, we have determined the phase diagram and phase transition 
behavior of the extended model. Our results are fully consistent with scal- 
ing theory (for the case when a "marginal" field is present) despite the 
unusual nature of the transition for h = 0. In particular, we find for the 
renormalization group eigenvalues yh — Ut — d, and for the sub-leading 
eigenvalues z t — x and — 0. We also examine a completely solvable 
model with very similar thermodynamics, but for which all phase transi- 
tions are first-order. 
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APPENDIX: BOUNDS FOR %±l 

First we introduce some notation (following [2|). 

We use rffl := for the fractions (called Farey fractions), where n 

is the order of the Farey fraction in level N. Level N = consists of the 
two fractions {°,|}. Succeeding levels are generated by keeping all the 
fractions from level N in level N + 1, and including new fractions. The 
new fractions at level N + 1 are defined via g?^™\ := (%' + e?^ +1 ^ and 

n N+i := n 7V^ + n N +1 ^ so tnat 
N = {?,!} 

N = 2 |S 1 I 2 11 etc 

Note that n = 1, . . . , 2 N + 1. When the Farey fractions are defined using 
matrices (spin states) A and B, the level N + 1 is the number of matrices 
in the chains starting with matrix A and hence the length of the spin chain 

Q . 

Using this notation we can write the partition function @ restricted to 
chains starting with A 

Z${p,h) = —r* r —. — , /3 em. (A.l) 

Note that the partition function © is the sum of Z^(f3, h) and Z^(P, h), 
where the Z^((3,h) is the partition function for chains starting with the 
matrix B. First we find bounds for Z^(f3, h) and then prove a lemma which 
lets us apply the bounds for Z^((3, h) to Z^(/3, h) also. 

Now, when we go from level N to level N + 1 we double the number 
of the terms in the partition function. Note that for chains starting with 
the matrix A one half of the terms come from matrix products of the form 
AMm-iA and the others from products AM^-iB. It is easy to check 
that the corresponding traces for given n £ {1, . . . , 2 N } are <%+7^ + n jv+i 
and <i^™\ + H^T^ 1 "*, respectively. These traces are multiplied by an h 
dependent factor e^ /3h ^ 2 ^-' i = 1 t7i ^ N ^ 1 ) which is simply e@ h raised to the 
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power (#A — #B), the number of matices A minus the number of matices B 
in the particular chain. For the terms from products of the form AM^-iA, 
it follows on using the definition of the Farey fractions that 

p -fJh(2-£t\ 1 ^- N - 1 ) p -Ph(2Z" = i^- N )+f }h e -/"»(2£!Li<7i--N) 

< P P\ h \ 



Ofc 1 ' + n%l\y (4 l) + ^ + ri" +1 V " (d£> + n^Y 
and, similarly, for AM^-iB 

e -l3h(2Z" = \ 1 <Ti-N-l) e -[3h(2Z? =1 ^-N)-f3h e -0h(2Y;? =1 <Ti-N) 

(ST+^SfV = (4" ) + 4" +1) +4 +1) )^ " (4 } +4 l+1 V e 

For the lower bound we just need the AMn-iA terms 

e -/3/i(2E 4 =i^-JV)+/Jh e -Ph(2Y:ti^-N) £ -/3\h\ 



where we used the fact that nffl < dffi. Thus we get for Z^ +1 (j3,h) = 

z^r iA (&h)+z™r iB (p,h) 

2-ee-^z£(f3,h) < Z$ +1 {(3,h) < 2e?Wz£{(3,h) 

for any > and h E R. 

Finally, we prove a lemma which allows us to bound Z^ (f3,h). Consider a 

(2 x 2) matrix M = ( mi ™ 2 ) and define the operator ~ via M := ( mi ms ) . 

Then we have the following result. 



Lemma.l Let M = AZiZ 2 . . . Z N , where Z { e {A, B}, with A 



:;?) 

and B = Then M — BZ\Z 2 ■ ■ -Zjq, i.e. the ~ operator exchanges 

A and B. 



Proof. We will use mathematical induction. It is easy to see that A = B 
and B = A. From matrix multiplication follows BM = ( m2 „j" 14 mi +™ 3 ) 

andAM=( m 2 ™ 2 V 

Clearly the ~ operation is a 1-to-l map of the set of all chains AMn 
onto BMm- Furthermore, the magnetic field term in the energy of each 
chain changes sign under this operation, so that the bounds just obtained 
for Zjf(f3, h) may be applied to Z§(/3, h). Therefore 

Zn 
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Note that the proof is easily adapted to the KSC model. 
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